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ABSTRACT 

We present spectroscopic observations in the rest-frame optical and near- to mid-infrared wavelengths of four 
gravitationally lensed infrared (1R) luminous star-forming galaxies at redshift 1 < z < 3 from the LUCIFER 
instrument on the Large Binocular Telescope and the Infrared Spectrograph on Spitzer. The sample was selected to 
represent pure, actively star-forming systems, absent of active galactic nuclei. The large lensing magnifications result 
in high signal-to-noise spectra that can probe faint IR recombination lines, including Paa and Bra at high redshifts. 
The sample was augmented by three lensed galaxies with similar suites of unpublished data and observations from 
the literature, resulting in the final sample of seven galaxies. We use the IR recombination lines in conjunction 
with Ha observations to probe the extinction, A„, of these systems, as well as testing star formation rate (SFR) 
indicators against the SFR measured by fitting spectral energy distributions to far-IR photometry. Our galaxies 
occupy a range of A v from ~0 to 5.9 mag, larger than previously known for a similar range of IR luminosities at 
these redshifts. Thus, estimates of SFR even at z ~ 2 must take careful count of extinction in the most IR luminous 
galaxies. We also measure extinction by comparing SFR estimates from optical emission lines with those from far- 
IR measurements. The comparison of results from these two independent methods indicates a large variety of dust 
distribution scenarios at 1 < z < 3. Without correcting for dust extinction, the Ha SFR indicator underestimates 
the SFR; the size of the necessary correction depends on the IR luminosity and dust distribution scenario. Individual 
SFR estimates based on the 6.2 pm polycyclic aromatic hydrocarbon emission line luminosity do not show a 
systematic discrepancy with extinction, although a considerable, ~ 0.2 dex, scatter is observed. 
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1. INTRODUCTION 

The evolution of the star formation rate (SFR) of galaxies is 
a central topic to the study of galaxy evolution. It is generally 
agreed that the SFR density of the universe has declined by 
an order of magnitude since z ~ 1 to the present (e.g., Le 
Floc’h et al. 2005; Rujopakarn et al. 2010; Magnelli et al. 201 1). 
The exact epoch of the peak of the SFR history is not known 
precisely, although it appears to be constrained to be within 
1 < z < 3 (e.g., Perez-Gonzalez et al. 2005; Reddy et al. 2008; 
Rodighiero et al. 2010; Magnelli et al. 201 1). 

At this redshift range, the primary SFR indicators are based on 
infrared (IR), optical, and extinction-corrected ultraviolet (UV) 
observations; the resulting SFR estimates commonly disagree 
with each other by more than a factor of two (e.g., Reddy 
et al. 2008). The majority of star formation at these redshifts 
is known to occur in optically extincted star-forming regions in 
IR luminous galaxies (e.g., Le Floc’h et al. 2005; Dole et al. 
2006; Berta et al. 2011), and the uncertainties resulting from 
the extinction undermine our ability to understand the SFR 


* The LBT is an international collaboration among institutions in the United 
States, Italy and Germany. LBT Corporation partners are: The University of 
Arizona on behalf of the Arizona university system; Istituto Nazionale di 
Astrofisica, Italy; LBT Beteiligungsgesellschaft, Germany, representing the 
Max-Planck Society, the Astrophysical Institute Potsdam, and Heidelberg 
University; The Ohio State University, and The Research Corporation, on 
behalf of The University of Notre Dame, University of Minnesota, and 
University of Virginia. 


history during this era. A better understanding of the SFR history 
will have important implications on the cosmic stellar mass 
buildup and metal production. For example, there is a significant 
discrepancy between the expected metal abundance derived 
from the SFR history and the observed abundance (Bouche 
et al. 2007), with some studies indicating the difference to be 
nearly an order of magnitude, and as a result placing the peak 
of the SFR history as far back as z ~ 4 (Kobayashi et al. 
2007). 

The situation necessitates more efforts into exploring the 
nature of optical extinction in high-redshift galaxies using 
unbiased measures as well as studying the effect of extinction 
on various SFR indicators. 

Extinction measurements based on optical emission lines 
could be biased in highly obscured environments because 
the indicator only probes the outer layer of the star-forming 
regions, where the extinction is relatively low. This effect 
is observed locally (e.g., Alonso-Herrero et al. 2006), but 
only because observations are available for longer wavelength 
IR recombination lines (e.g., Paa and Bra) that are less affected 
by extinction. In this work, we extend this technique out 
to z — 3 by comparing the strength of the Ha line with 
those of Paa and Bra. Since the latter are in wavelength 
regions with 10-20 times less extinction than Ha, they provide 
a measurement of extinction through the entire star-forming 
region. With an unbiased estimate of extinction, we can address 
its effects on SFR indicators in this critical redshift range. 
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In this paper, we study a wide range of star formation diag- 
nostics in seven gravitationally lensed star-forming galaxies at 
1 < z < 3. We observed four of these galaxies spectroscopically 
with the LUCIFER 1 near-IR imaging spectrograph on the Large 
Binocular Telescope (LBT) to measure the Ha line flux in the 
near-IR, and with the Infrared Spectrograph (IRS, 5-38 pm) on 
the Spitzer Space Telescope to observe the wavelength regions 
covering Paa and / or Bra lines as well as aromatic emission lines 
and emission complexes (commonly attributed to, and hereafter, 
polycyclic aromatic hydrocarbons or PAHs). The sample com- 
prises A2218b, A2667a, A2218a, and A1835a at redshift 1.03, 
1.03, 2.52, and 2.57, respectively (magnification 6x-27x). The 
IRS observations at longer wavelengths are further used to com- 
pare SFR estimates from the 6.2 pm PAH feature. They were 
carried out under Spitzer program ID 82, 30775, and 50586 (PI: 
G. Rieke); and 40443 and 50372 (PI: C. Papovich). Our sample is 
augmented by three galaxies with similar suites of observations 
from unpublished data and the literature: SDSSJ 12060 1+5 142 
(hereafter the Clone), the Lyman break galaxy LBG MS 1512- 
cB58 (hereafter cB58), and the 8 O’clock arc, at redshifts 2.00, 
2.73, and 2.73, respectively. Our galaxies were also observed 
with the MIPS instrument (Rieke et al. 2004) on Spitzer at 24 
and 70 pm by Rigby et al. (2008), allowing us to combine the 
MIPS data with far-IR photometry from the literature to estimate 
SFR via L(TIR) obtained by fitting the spectral energy distribu- 
tion (SED) across the peak of the far-IR emission. In addition, 
we used the 24 pm monochromatic indicator (Rujopakarn et al. 
2012) to estimate L(TIR). 

This paper is organized as follows. We discuss the LBT and 
Spitzer observations and data reduction in Section 2; extinction, 
metallicity, L(TIR), and SFR measurements in Section 3; 
compare SFR indicators in Section 4; and make concluding 
remarks in Section 5. We assume a ACDM cosmology with 
Q. m — 0.3, Ll\ = 0.7, and Hq = 70kms _1 Mpc - 1 throughout 
this paper. 

2. OBSERVATIONS AND DATA REDUCTION 

Here we describe the selection of our sample (Section 2.1); the 
spectroscopic observations with LBT /LUCIFER along with the 
near-IR data reduction (Sections 2.2 and 2.3) and with Spitzer 
(Section 2.4); data compilation from the literature (Section 2.5); 
and additional photometric observations with Spitzer MIPS and 
the respective data reduction procedure (Section 2.6) that results 
in the final data set for our sample of galaxies. 

2.1. Sample Selection 

To measure extinction and study its implication on the 
SFR indicators, we need a sample of star-forming galaxies 
that we understand very well. Gravitationally lensed galax- 
ies are outstanding candidates because their amplified fluxes 
and images enable high signal-to-noise (S/N) spectroscopy. 
The candidates for our new observations were drawn from the 
pool of objects studied by Rigby et al. (2008). Briefly, the Rigby 
et al. (2008) selection requires that objects have (1) Spitzer 
MIPS 24 pm flux above 0.4 mJy to be observed efficiently with 
IRS, (2) morphologies in the optical Hubble Space Telescope 
imaging that exclude members of the lensing cluster, (3) spec- 
troscopic or probable photometric redshift above 1, and (4) 
lensing magnification above 3 x . From the Rigby et al. (2008) 
sample, we further require that candidates are absent of active 
galactic nucleus (AGN) activity based on X-ray luminosity be- 
low 10 42 erg s -1 cuT 2 and the lack of mid-IR power-law SEDs, 


Table 1 

Sample of High-redshift Strongly Lensed Star-forming Galaxies 


Source 

z 

R.A. 

(J2000) 

Decl. 

(J2000) 

Lensing 

Magnification 

Reference 

A2218b 

1.034 

16 35 55.16 

+66 11 50.8 

6.1 

1,2 

A2667a 

1.035 

23 51 40.00 

-26 04 52.0 

17 

1,2 

The Clone 

2.003 

12 06 01.80 

+51 42 30.7 

27 ± 1 

3,5 

A2218a 

2.520 

16 35 54.18 

+66 12 24.8 

22 ±2 

1,2 

A1835a 

2.566 

14 01 04.96 

+02 52 24.8 

3.5 ±0.5 

1,2,6 

cB58 

2.729 

15 14 22.29 

+36 36 25.7 

~30 

4,7 

8 O’clock 

2.731 

00 22 40.97 

+14 31 14.0 

~8 

8,9 


References. (1) redshifts from our LBT Hof spectroscopy; (2) Rigby et al. 
2008; (3) Hainline et al. 2009; (4) Teplitz et al. 2000; (5) Lin et al. 2009; 
(6) Smail et al. 2005; (7) Seitz et al. 1998; (8) Allam et al. 2007; (9) Finkelstein 
et al. 2009. 

indicating that their 24 pm fluxes are not dominated by AGN 
emission. Apart from the objects selected from the Rigby et al. 
(2008) sample, we include three additional objects that match 
the selection criteria: cB58, the Clone, and the 8 O’clock arc. 
We will discuss these objects in Section 2.5. 

The galaxies in our sample are located behind galaxy clusters 
that are well modeled for mass distribution. The mass model of 
the A2218 cluster (Kneib et al. 1996) indicates a magnification 
value for A22 1 8a of 22 ± 2 (Kneib et al. 2004). For A22 1 8b and 
A2667a, magnification estimates of 6.1 and 17, respectively, 
are modeled by Rigby et al. (2008). For A1835a, we adopt the 
magnification estimate from Smail et al. (2005), 3.5 ± 0.5. The 
sample is listed in Table 1. 

2.2. LBT /LUCIFER Near-infrared Spectroscopy 

Our near-IR spectroscopy was obtained in three campaigns 
during 2010 October and 2011 May using the LUCIFER 1 
instrument on the LBT (Mandel et al. 2007). The common 
aspects of all these observations are that we use the N 1 .8 camera 
with plate scale of O'.' 25 pixel -1 and with the 210 lines mm -1 
high-resolution grating, resulting in spectral resolution of 5.0 A. 
A-stars were observed for spectrophotometric calibration for 
each target (details on spectral type below) at similar airmass 
as the science exposure, either immediately before or after the 
science exposures. Positions of the LUCIFER slit for each of 
our observations are shown in Figure 1 . 

A2667a was observed on 2010 November 3 in the long-slit 
spectroscopy mode using a slit of 1/5 by 3.9 arcmin under 
variable seeing of 1/2-2/5. The slit was rotated to P.A. = 320° 
to place it along the length of the lensed arc. The total integration 
time was 12 x 300 s with the telescope nodded 50" along the slit. 
The unusually long nodding was required by the ~25" length 
of the arc. For calibration, we observed HD 223466, an A3V 
calibration star with V — 6.42 mag, through the same setup using 
a total integration time of 2 x 20 s in two nodding positions. 

A2218b was observed on 2011 March 10 in the long-slit 
spectroscopy mode using a slit of I/O by 3.9 arcmin under 
0/ 4-0/5 seeing. The slit was rotated to P.A. = 42/5 to place it 
along the lensed arc, which is ~20" long. The arc is 1/9 wide in 
the widest region; thus, we took the spectra at the central position 
and another position shifted by 0/9, with the 0/ 1 overlapping 
allowing for possible pointing errors to prevent a gap in the slit 
mapping. We will refer to these positions as the “Central” and 
the “Central + 0/9” positions, respectively. The total integration 
time was 6 x 300 s at the Central position and 5 x 300 s at the 
Central + 0/9 position, totalling 55 minutes, with the telescope 
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Figure 1 . Images of the lensed galaxies in our sample from the Hubble Space Telescope Advanced Camera for Survey (ACS) with the positions of our LBT /LUCIFER 
slit for NIR spectroscopy shown in red. Clockwise from top left are A2218b, A2218a, A1835a, and A2667a. The ACS filter for A2218a and A2218b is F625W; A1835 
and A2667 images were taken in F850LR Each image is 37" wide (identical scale in all four panels), north is up, east is left. For A2218a and A2667a, each lensed 
component is labeled following Kneib et al. (2004) and Sand et al. (2005), respectively. 

(A color version of this figure is available in the online journal.) 


nodded 40" along the slit for both. An A0V calibration star, 
HD 145454, with V = 5.43 mag was observed using a total 
integration time of 2 x 20 s in two nodding positions. 

A2218a was observed on 2011 March 11 in the long-slit 
spectroscopy mode using a slit of l'.'O by 3.9 arcmin under 
O'.'9-l'.'l seeing. The slit was rotated to P.A. = 18°9 to place it 
along the lensed arc. The total integration time was 10 x 300 s 
with the telescope nodded 20" along the slit. We observed the 
same calibration star as for A22 1 8b. 

A1835a was observed on 2011 May 7 in the long-slit 
spectroscopy mode using a slit of 1"5 by 3.9 arcmin under 
l'.'O-l'.'l seeing. The slit was rotated to P.A. = 15°0 to place it 
along the length of A1835a. The total integration time was 12 
x 300 s with the telescope nodded 20" along the slit. An A2V 
calibration star, HD 122365, with V — 5.98 mag was observed 
using a total integration time of 2 x 20 s in two nodding positions. 

We took Ar lamp exposures for wavelength calibration along 
with dark and flat exposures during daytime. 

2.3. LBT/ LUCIFER Data Reduction 

Our near-IR spectral reduction has four steps: sky subtraction, 
spectral extraction, flux calibration, and line flux measurements. 
We first use a modified version of an IDL long- slit reduction 
package written by G. D. Becker for NIRSPEC (Becker et al. 
2009) to perform sky subtraction following the prescription of 
Bian et al. (2010). The sky-subtracted two-dimensional spectra 
are corrected for the two-dimensional dispersion distortions, 
wavelength calibrated, and extracted using IRAF. We note that 
the flexure compensation mechanism of LUCIFER1 was not 


available during our campaigns, but clean sky subtraction was 
achieved with the reduction procedure. We combine spectra for 
each object by averaging values within each wavelength element 
with sigma clipping at 2.5cr and then flux calibrate using the 
xtellcor_general tool (Vacca et al. 2003) and the A-star 
calibrators described in the previous section. Line fluxes for our 
objects are measured by integrating the line and subtracting 
the background estimated from the mean continuum in the 
spectral range. We estimate the line flux uncertainties by a Monte 
Carlo simulation done by repeatedly (n = 10 4 ) simulating the 
object’s spectrum using random values drawn from the Gaussian 
distribution centered at the observed flux and a equal to the 
observed uncertainties and repeat the line measurement. We take 
the resulting la distribution of the simulated line flux values to 
be the uncertainties for the line flux. 

2.4. Spitzer IRS Mid-infrared Spectroscopy and Data Reduction 

Unpublished Spitzer IRS spectroscopy in this work is from 
(1) Spitzer program 50586 (PI: G. Rieke) to take deep (in- 
tegration times ~15ks) first-order short-low spectra (SL1; 
7.4-14.5/xm) in the Bra wavelength regions of A2667a and 
A2218b; and (2) Spitzer program 40443 (PI: C. Papovich) to 
take deep (integration time 7-14 ks) second-order short-low 
spectra (SL2; 5. 2-7. 7 /xm) in the Paa wavelength regions of 
A2218a and A1835a. For each of the four objects, the IRS 
first-order long-low spectra (LL1; 19.5-38.0 /xm) were from 
GTO programs 82 (A2218a and A 1835a) and 30775 (A2218b 
and A2667a), PI: G. Rieke. The IRS second-order long-low 
spectra (LL2; 14.0-21.3 /xm) are available for A2218b and 
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Figure 2. Spitzer / IRS mid-IR and LBT/LUCIFER near-IR spectra of A2218b. Panel (a) is the magnified region of IRS SL1 spectrum at the Brel' wavelength. The 
best Gaussian fit to the Bra emission line is the thick black line. The color-coded Gaussians at the Bra line wavelength illustrate the expected Bra line flux for each 
extinction scenario from A„ of 0 to 8.0 mag given the observed Ha flux. Panel (b) shows the Ha emission line from LBT/LUCIFER. The light and dark gray spectra 
are extracted from the two parallel slit positions shown in Figure 1 . A double-peak line profile is clearly visible. The continuum level for each aperture is shown 
as dashed lines. Panel (c) presents the overall IRS spectra with the green, orange, and red lines representing spectrum from the IRSSL1, LL2, and LL1 modules, 
respectively. 

(A color version of this figure is available in the online journal.) 


Table 2 

Near- and Mid-infrared Spectroscopic Observations 


Source 

SL2 

IRS Exposure Time (ks) 
SL1 LL2 

LL1 

IRS Program ID 

Near-IR Spectroscopy 

A2218b 


15.1 

3.6 a 

3.6 a 

30775, 50586 

This work 

A2667a 


16.1 

1.9 a 

1.9 a 

30775, 50586 

This work 

The Clone 

8.6 


21.6 b 

20.9 b 

40430, 50372 

Hainline et al. (2009) 

A2218a 

14.4 C 



3.6 a 

82, 40443 

This work 

A1835a 

7.7 



3.6 a 

82, 40443 

This work 

cB58 

6.8 d 

6.8 d 

14.6 C 

34. l c 

30832 

Teplitz et al. (2000) 

8 O’clock 

14.4 



22.1 

40443 

Finkelstein et al. (2009) 


Note. Spitzer IRS spectra previously published in (a) Rigby et al. (2008), (b) Fadely et al. (2010), (c) Papovich et al. 
(2009), and (d) Siana et al. (2008). 


A2667a from program 30775, and A1835a from program 82. 
These data are published in Rigby et al. (2008). The details of 
IRS observations for our sample are summarized in Table 2. 

We have reduced all archival spectroscopic data for objects 
in our sample using the most updated version of the processing 
software from Level- 1 data (BCD) for IRS by following the 
prescription of Teplitz et al. (2007). Briefly, first we use IRAF 
to fit the background slope with time for each row and subtract 
that fitted value row by row to remove latent charge. Second, we 
interpolate over bad and/or hot (“rogue”) pixels using the IDL 
routine IRSCLEAN (version 2.1) provided by the SSC based on 


the known hot pixel mask for the corresponding campaign, plus 
manual identification. Third, we subtract sky using a sky image 
constructed from a median of other map positions. Then we co- 
add each map position into a two-dimensional (2D) spectrum 
and use the SPICE software (version 2.4) provided by the SSC 
to extract a one-dimensional (ID) spectrum using the optimal 
extraction template for point sources for each map position. The 
ID spectra are combined with sigma clipping at 2.5er into the 
final ID spectra shown in Figures 2 and 3. 

We use the PAHFIT package (Smith et al. 2007) to measure 
PAH feature fluxes. PAHFIT uses / 2 minimization to simulta- 
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Figure 3. Same as Figure 2 but for A2667a, A2218a, and A1835a, from top to bottom, respectively. The IRS SL2 module spectra (not available for A2218b in the 
previous figure), which cover Paa lines for A2218a and A1835a, are shown in blue. There is no SL1 and LL2 coverage for A2218a and A1835a. 

(A color version of this figure is available in the online journal.) 


neously fit the PAH and nebular emission lines, stellar and dust 
continua, as well as the silicate absorption feature. We report the 
PAH flux measurements in Table 3. To define the dust and stellar 
continuum for the measurement of Paa and/or Bra line fluxes, 
we iteratively fit the IRS spectra with PAHFIT and interpolate 
over points with fit residuals greater than 3cr. This process is 
repeated until there are no remaining outliers (usually within 
the third iteration) and the final fit, which is effectively a noise- 
free template representing the object’s spectrum, is adopted as a 
continuum estimate. The actual line flux measurements for both 
the PAH emissions and recombination lines are done on the orig- 
inal spectra, i.e., only the continua for IR recombination line 
measurements are determined with this iterative fitting proce- 
dure. This method of continuum estimation allows construction 
of a mid-IR spectral template that best matches each galaxy 
in a non-parametric manner. In other words, it provides con- 
tinuum estimates in the Paa and Bra wavelength regions that 


are not only constrained by photometric observations (e.g., from 
Spitzer/IRAC), but also by PAH-region information. These con- 
tinuum estimates are shown for each object in Figures 2-4. 

Bra and Paa lines are measured by fitting Gaussian profiles 
with widths fixed at the resolution of the IRS module covering 
the line, the central wavelengths fixed at the expected line 
wavelengths based on the optical spectroscopic redshift, and 
the Gaussian base given by the continuum under the expected 
wavelength range of the line; the line peak is the only free 
parameter to fit. The line flux uncertainties are estimated by 
a procedure similar to that of the Ha line flux uncertainties: 
(1) a series of spectra (n = 10 4 ) were generated randomly 
based on the actual spectra and the a uncertainties, (2) line 
flux measurements are carried out on them, and (3) we take the 
ler value of the simulated flux distribution to be the uncertainties 
for the line flux. If the formal fit to the line yields negative line 
flux due to non-detection (occurring in two cases, cB58 and the 
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Rest Wavelength (/x m) 

1.8 1.9 2.0 2.1 2 4 6 8 10 12 



Observed Wavelength (/xm) 

Figure 4. Our reduction of the Spitzer IRS spectra for the Clone (SDSSJ 1206+5 142) and cB58 (cf. Fadely et al. 2010 and Siana et al. 2008, respectively). The color 
coding for IRS modules the line fluxes for each A v scenario is the same as Figures 2 and 3. We note that the apparent spectral feature at ~5.9 /xm (observed frame) 
in the SL2 spectrum of the Clone cannot be interpreted as a Paa line at the redshift z = 2.00, which has been measured from multiple optical lines (see discussion in 
Section 3.1). 

(A color version of this figure is available in the online journal.) 


Clone), we adopt the 1 a value as an upper limit to the line flux. 
The measurements are reported in Table 4. 

2.5. Archival and Literature Data 

We augmented our sample with three additional objects: the 
Clone, cB58, and the 8 O’clock arc. These objects have suites 
of data similar to our sample: rest-frame optical spectroscopy 
covering Ha and deep IRS spectroscopy covering Paa (Table 1). 

cB58 (Yee et al. 1996), at redshift z = 2.729, is the first 
lensed LBG known with high magnification (~30x); it has 
been studied extensively from UV to millimeter wavelengths 
(Siana et al. 2008, and references therein). Rest-frame optical 
spectroscopy was obtained by Teplitz et al. (2000). Siana et al. 


(2008) present a full suite of IRS observation (SL2, SL1, LL2, 
and LL1) for the galaxy. We re-reduced these IRS data using our 
procedure above and found the result to be consistent with the 
original reduction published in Siana et al. (2008); details are 
given in Section 3.1. The lensing magnification of cB58, ~30x, 
was from Seitz et al. (1998). 

The Clone, at redshift z = 2.0026, was discovered 
from SDSS imaging (Lin et al. 2009). Its rest-frame op- 
tical spectroscopy was obtained by Hainline et al. (2009). 
IRS observations in the LL2 and LL 1 bands were obtained by 
Fadely et al. (2010), who concluded that the Clone’s IR emission 
is dominated by a starburst based on the strength of the 6.2 /xm 
PAH feature, despite the object showing a very strong [Siv] 
emission line. The Clone’s deep IRS SL2 observation is a part 
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Source 

/(PAH 6 .2) 

/(PAH77 Complex) 

/(PAH 86 ) 

/(PAH n 3 Complex) 

/(PAH 12 . 6 ) 

A2218b 

15.9 ± 0.2 

72.2 ± 1.8 

13.2 ± 1.0 

17.3 ± 0.7 

12.7 ± 0.9 

A2667a 

13.0 ±0.2 

82.1 ± 3.0 

19.9 ± 2.2 

21.6 ± 1.3 

11.7 ± 1.6 

The Clone 

7.5 ± 0.3 

15.9 ±0.4 

4.2 ± 0.3 

4.5 ± 0.7 


A22 18a 

16.0 ± 0.8 

54.4 ± 1.3 

10.3 ±0.8 



A1835a 

7.5 ± 1.4 

26.1 ± 1.6 

7.5 ± 1.9 



cB58 

2.9 ± 0.2 

9.7 ±0.5 

2.5 ± 0.4 




Notes. PAH emission lines and emission complex fluxes measured by PAHFIT in 10 -15 ergs _1 cm - ". These 
values are not corrected for lensing magnification. Dots indicate that the feature is outside spectral coverage. The 
7.7 (tm and 11.3 // in complexes are the sum of fluxes at wavelength 7. 3-7. 9 /xm and 1 1.2-11.4 /xm, respectively. 
The LL2 spectrum of the 8 O’clock arc does not have sufficient S/N to measure PAH line fluxes consistently. 


Table 4 

Measured Fluxes 


Source 

fl 24 Mm) 
(1) 

fOOum) 

(2) 

/(Ha) 

(3) 

/([Nil] 6583) 
(4) 

/( Paa) 
(5) 

/(Bra) 

(6) 

/[Ne 11 12.8] 
(7) 

/[Ne hi 15.5] 
( 8) 

A2218b 

1.67 

7.4 ± 1.5 

21.2 ±0.6 

5.2 ± 0.5 


4.3 ± 0.9 

1.8 ± 0.2 

0.9 ± 0.2 

A2667a 

4.52 

19.4 ±3.9 

78.7 ± 1.5 

31.1 ± 3.2 


0.4 ± 1.0 

2.8 ± 0.4 

0.9 ± 0.3 

The Clone 

0.88 

8.8 ± 1.8 

20.2 ± 0.5 a 

3.9 ± 0.2 a 

<1.3 




A2218a 

1.16 

<7 

6.2 ±0.3 

1.5 ±0.3 

5.7 ± 1.2 




A1835a 

0.99 

<13 

1.8 ±0.3 

1.0 ±0.3 

7.9 ± 1.5 




cB58 

0.24 b 

1.7 ± 1.0 b 

12.6 ± 0.4 C 

1.1 ±0.3 C 

<2.3 

<1.3 



8 O’clock 

0.57 


17.6 ± 0.5 d 

4.8 ± 0.4 d 

11.5 ±2.2 





Notes. Column 1 : 24 gm flux in mJy./(24 //m ) of Abell galaxies except A2667a are from Rigby et al. (2008), who noted that DAOPHOT errors 
of the flux are overly optimistic and thus not reported. A nominal uncertainty of 0.1 mJy is adopted for SED fitting for L(TIR). Rigby et al. 
(2008) mistakenly reported an/(24 // m ) for only one component of the A2667 arc. A correct value is reported here. Column 2: 70 // m flux in 
mJy. Columns 3 and 4: optical line fluxes in 10 16 erg s -1 cm -2 . Columns 5 and 6: IR recombination lines fluxes in 10 16 erg s cm -2 . The 
formal fit to the Bra line of cB58 yields a flux value of (—0.4 ± 1.3) x 10 16 erg s' cm -2 and hence a 1 a limit is reported in the table. Also for 
the Pact line, the formal fit yields (—0.2 ± 2.3) x 10’ 16 erg s -1 cm -2 and hence the lcr limit is reported. Likewise, a Iff limit is reported for the 
Paa line flux of the Clone, where a formal fit yields (0.3 ± 1.3) x 10 -16 . Columns 7 and 8: line fluxes of [Ne n] 12.8 /cm and [Ne hi] 15.5 (im 
in if) 15 erg s' 1 cm -2 measured by PAHFIT. Dots indicate no spectral coverage. These values are not corrected for lensing magnification. Flux 
references: (a) Hainline et al. (2009); (b) Siana et al. (2008); (c) Teplitz et al. (2000); (d) Finkelstein et al. (2009). 


of Spitzer program 50372 (PI: C. Papovich) with an integration 
time of 8.6 ks. As with cB58, we re-reduced all IRS data using 
our procedure and software to ensure consistency. We present 
our reduction of these archival data in Figure 4. The lensing 
magnification for the Clone, 27 ± 1 x , was reported by Lin et al. 
(2009). 

The 8 O’clock arc at redshift z = 2.7308 was discovered 
serendipitously from SDSS imaging by Allam et al. (2007). 
Finkelstein et al. (2009) obtained the rest-frame spectroscopy 
of components “A2” and “A3” in the nomenclature of Allam 
et al. (2007). The arc was observed as a part of Spitzer program 
40443 (PI: C. Papovich) with an SL2 integration time of 14.4 ks. 
The slit for the SL2 observation is centered on the “A2” 
component and covered significant parts of both the “Al” and 
“A3” components such that any bias due to spatial variation of 
recombination lines along the lensed arc will be negligible. The 
PAH spectrum of the 8 O’clock arc shows 6.2 /xm and 7.7 /xm 
emission features, yet at low S/N due to contamination in the 
IRS slit from a nearby IR-bright dust-obscured galaxy, which 
we will discuss in the next section. We adopt the magnification 
estimate of 8x for the arc (Finkelstein et al. 2009). 

2.6. Additional Spitzer MIPS Photometry and Data Reduction 

The Spitzer MIPS 24 and 70 /xm photometry for objects 
behind the Abell clusters are taken from Rigby et al. (2008) with 
exception of A2667a, whose 70 /xm photometry is reported as 
an upper limit in Rigby et al. (2008) but which was re-observed 
with longer integration time in Spitzer program ID 50586 (PI: 


G. Rieke). For the Clone, we downloaded the 24 and 70 /xm 
data observed as part of programs 40430 (PI: S. Allam) and 
50372 (PI: C. Papovich). The Spitzer MIPS observations for 
the 8 O’clock arc were obtained under program 40443 (PI: C. 
Papovich). 

We measured the 24 and 70 /xm fluxes from the Level- 
2 (PBCD) data using the Starfinder IDL routine (Diolaiti 
et al. 2000) to perform PSF photometry. The PSF model was 
generated using the STinyTim routine provided by the SSC and 
smoothed according to a prescription given by Engelbracht et al. 
(2007) to better match the observed PSF. We require o = 3 
and the PSF correlation value of 0.75 for our 24 and 70 /xm 
flux measurements. The 8 O’clock arc has a 24 /xm bright dust- 
obscured galaxy ~6" 5 away that partially blends with the arc, but 
we were able to deblend the flux by performing PSF photometry 
to subtract the intervening source first and repeat the photometry 
to measure the flux of the arc. This contaminating galaxy is 
highly extincted, such that it is undetected in the SDSS imaging 
and was discovered only in the 24 /xm observation. Our flux 
measurements along with fluxes from the literature are reported 
in Table 4. 

2. 7. Augmented Sample Summary 

Our augmented sample consists of the four galaxies we ob- 
served plus the three from the literature. In addition to the mea- 
surements of Ha and of at least one of the IR recombination 
lines (Paa and/or Bra), all of them have good 24 /xm pho- 
tometry, far-IR and sub-millimeter photometry as well as 
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IRS spectroscopy covering, at least, the 6-8 pm wavelength re- 
gion. The IRS coverage goes out to nearly 17 /an (rest frame) in 
the two lowest redshift objects (A2667a and A2218b, z ~ 1.0). 

3. RESULTS 

In this section, we describe the near- and mid-IR spectra and 
the derived physical quantities, including metallicity, L(TIR), 
and SFR from optical and IR tracers, and extinction measure- 
ments. 

3.1. Near- and Mid-IR Spectroscopy 

The final ID Ha spectra from our LBT/LUCIFER1 observa- 
tions are shown in Figures 2 and 3 along with the IRS spectra. 
The 6.2 /im PAH emission line, 7.7 /Tin complex, and 8.6 pt m 
complex are clearly detected in all galaxies. 

The spectra from both slit positions (Section 2.2) of A2218b 
show highly asymmetric line profiles. The profile at the cen- 
tral position is double-peaked in both Ha and [N ii] while the 
Central + O'.' 9 spectrum shows a single-sided profile with 
stronger emission at the red side of the line. The combined 
line profile from both slit positions is well fitted with a dou- 
ble Gaussian, whose deconvolved spectral width (corrected 
for instrumental resolution) indicates a velocity dispersion of 
163 ± 3kms _1 . Traces of the arc northeast of the main com- 
ponent hinted in the ACS image are not detected at Ha. The 
IRS spectrum for A2218b shows a strong [Nem] 15.56 /tm 
line as well as the [Nen] line at 12.81 pm, apart from other 
PAH features. The Bra line is detected at 4.8a with line flux of 
(4.3 ± 0.9) x 10” 16 ergs _1 cm” 2 . 

The spectra of A2667a were extracted from two apertures 
in the slit. Both of these positions are very bright in Ha. 
Referring to Figure 1, these apertures are at (1) the main 
arc north-northeastern of the brightest cluster galaxy (BCG, 
identified as “Al” by Sand et al. 2005, and in Figure 1); and 
(2) the smaller arc east of the BCG (the “A2” arc in the Sand 
et al. study). Although it appears that a part of the smaller arc 
is outside the slit in the optical image, the Ha line is clearly 
detected (the light gray spectrum in Figure 3, top row, left 
column). Both apertures show the same profile for both Ha and 
[N ii] lines. No velocity shift between the two components 
is observed. The formal fit to the Bra line yields a flux of 
(0.4± 1.0)x 1 0” 16 , a non-detection. The galaxy has the warmest 
dust continuum among our samples, as seen in the rising SED 
toward the red end of the IRS spectrum in Figure 3 (cf. A22 1 8b’s 
continuum in Figure 2). However, the strong PAH emission is 
inconsistent with the system’s emission being dominated by 
an AGN in the mid-IR. Its log([N ii]/Ha) line ratio of —0.40 
alone (i.e., without the [Om] and H/J) does not unambiguously 
support or rule out the presence of a weak narrow-line AGN. 
But even with the presence of a weak, optically selected AGN, 
the EW of the 6.2 pm PAH emission could still serve as a good 
tracer of SFR because the suppression of PAH emission due to 
AGN is found to be limited to the nuclear region (Diamond- 
Stanic & Rieke 2010). 

The Ha line of A2218a also shows a double-peaked profile 
indicating a velocity dispersion of 188 ± 5kms _1 . The galaxy 
has been a subject of a previous near-IR spectroscopic study by 
Richard et al. (2011). Our slit orientation, however, is different 
(see their Figure 1); they excluded the (5 and y components 
(the former being a major component) in the nomenclature 
of Kneib et al. (2004), as shown in Figure 1, while our slit 
orientation covers all the main components (a, ft, and y). The 


Rest Wavelength (pm) 

1.80 1.85 1.90 1.95 2.00 



Figure 5. Spitzer IRS SL2 spectrum for the 8 O’clock arc, showing the Paa line. 
The color coding for IRS module and the line fluxes for each extinction scenario 
are the same as Figures 2-4. 

(A color version of this figure is available in the online journal.) 


Paa observation, which targets the ft component, detected Paa 
line at 4.8a with line flux of (5.7 ± 1.2) x 10” 16 erg s” 1 cm -2 . 
This measurement is consistent with the previously published 
reduction of the same data by Papovich et al. (2009), who 
measured the line flux to be (8.5 ± 1.4) x 10” 16 erg s” 1 cm -2 . 
We attribute the small difference in the two measurements 
to the different definition of the continuum: Papovich et al. 
(2009) defined their continuum with IRAC photometry whereas 
our reduction adopted the simultaneous multi-component fit of 
spectral features described in Section 2.4. 

A 1 835a shows a faint Ha line, too faint to measure asymmetry 
in the line profile, but a bright Paa line with a flux of (7.9 ± 
1.5) x 10” 16 erg s” 1 cm” 2 . This is the second brightest Paa 
detection in our sample, after the 8 O’clock arc, whose Paa flux 
is (11.5 ± 2.2) x 10” 16 erg s” 1 cm” 2 . TheLLl PAH spectrum of 
the 8 O’clock arc does not have sufficient S/N to use PAHFIT 
to construct a noise-free spectrum as a template for the Paa 
line flux measurement (Section 2.4) due to the contamination 
from the nearby IR-bright dust-obscured galaxy, thus instead 
the continuum of the arc for the measurement is defined by 
interpolating between the IRAC 5.0 and 8.0 pm photometric 
points. These IRAC observations yield fluxes of 74 ± 3 and 
64 ± 4 pSy for the 5.8 and 8.0 pm bands, respectively. We 
present the Paa spectrum of the 8 O’clock arc in Figure 5. 

For the two IRS spectra from the literature, the Clone and 
cB58, our reduction generally agrees with those published 
previously by Fadely et al. (2010) and Siana et al. (2008), 
respectively. Small differences are noted in the PAH line fluxes 
due to the way that the continuum for the fit is defined. PAHFIT 
uses multi-component dust features to represent the continuum, 
while a single-component power law is used by Fadely et al. 
and Siana et al. Their approach resulted in a similarly good fit 
to the data compared to ours, although the measured PAH flux 
could differ, as observed in the Clone’s PAH fluxes from the 
6.2-1 1.3 pm features where the Fadely et al. (2010) values are 
a factor of 2. 6-3. 7 larger than ours. Fadely et al. reported that 
their results are a factor of 1-8 higher than those of Brandi et al. 
(2006) in their re-reduction of the same data set. Brandi et al. use 
a continuum level defined on either side of the emission features 
and thus their PAH fluxes will be systematically smaller than 
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with either our PAHFIT measurement or the single-component 
continuum method. That is, we expect the measurements of the 
PAH line flux from the same data to yield flux values in the 
following increasing order: Brandi et al., PAHFIT, and Fadely 
et al./Siana et al., which is consistent with the differences 
we find. The formal Paa line fit for cB58 yields a flux of 
(—0.2 ± 2.3) x 10” 16 erg s” 1 cm -2 , from which we determined 
the Paa flux limit to be <2.3 x 10” 16 erg s” 1 cm” 2 , which 
agrees with the non-detection reported by Siana et al. (2008), 
<6 x 10” 16 erg s” 1 cm” 2 . The formal Bra line fit yields a line 
flux of (— 0.4 ± 1.3) x 10” 16 erg s” 1 cm” 2 , thus we report a limit 
of <1.3 x 10” 16 ergs _1 cm” 2 . 

The formal fit to Paa line flux of the Clone yields (0.3 ± 
1.3) x 10” 16 erg s” 1 cm” 2 , from which we report an upper limit 
similar to that of cB58 above. In the panel showing the Paa 
wavelength region for the Clone in Figure 4, we note a spectral 
feature at 5.9 /xm (1.96 /x m rest frame) that resembles a Paa line 
at j ~ 2.1. However, the spectroscopic redshift of 2.0026 was 
measured with multiple optical lines (Hainline et al. 2009) and 
thus precludes the possibility of the feature being interpreted 
as the Paa line. It is also unlikely to be an H2 emission line at 
1.96/xm given the absence of other H2 lines nearby (2.03 and 
2.12 /xm). The source of the apparent feature therefore remains 
unknown. 

3.2. Metallicity 

We estimate metallicity with the N2 index (Pettini & Pagel 
2004) using both linear and third-order fits (Table 6). The index 
is not affected by the optical extinction given the proximity of 
the lines. Both estimators yield consistent values, which agree 
well with those found by Erb et al. (2006) for massive galaxies at 
z > 2 selected via UV luminosity. The range of the metallicity 
found in our sample is also in good agreement with the sample 
of 28 lensed galaxies at 1.5 < z < 5 observed by Richard et al. 
(2011), 8.00 ^ Z ^ 8.94, with an average oxygen abundance 
of Z = 8.55, similar to the average value from the linear N2 
index in this work, Z = 8.56. That is, our sample has typical 
oxygen abundance for field and lensed galaxies at similar z. 

3.3. Infrared Luminosity and Star Formation Rate 

We can measure the SFR from L(TIR), which in turn can 
be estimated by three tracers: (1) far-IR SED fitting, (2) the 
monochromatic 24 /xm L(TIR) estimator given by Rujopakarn 
et al. (2012), and (3) PAH emission line luminosity. The SFR 
can also be estimated from the Ha luminosity as well as the 
combination of Ha and 24 /xm luminosity. We adopt the SFR 
derived from the SED-htted L(TIR) as a fiducial for comparison. 
SFRs are calculated from L(TIR) by the relationship given by 
Rujopakarn et al. (2012). All SFR estimates are converted from 
their respective IMF assumptions to the Kroupa (2002) IMF. The 
SFR results are presented in this section and their implications 
are discussed in Section 4. 

The far-IR SED fitting (e.g., 30 /< m to ~1 mm) is the 
most straightforward among the methods to measure L(TIR). 
However, far-IR photometry is often affected by confusion noise 
at longer wavelengths in crowded fields, thus limiting the 
applicability of far-IR SED fitting in deep cosmological surveys. 
The far-IR and submillimeter photometry in our sample is aided 
by magnification of the gravitational lenses and thus we can 
probe fluxes below the typical confusion limit of the instrument. 
The SED fitting for the L(TIR) for A2218b and A 1835a were 
carried out by Rigby et al. (2008) by combining the 24 and 


70 /xm data with submillimeter observations, including at least 
450 and 850 /xm for all their galaxies and additionally 1.3 mm 
for A1835a. The value for A2218a is measured using Herschel 
observations by Finkelstein et al. (201 1). The SED-fitted L(T1R) 
for the Clone and cB58 were from the Fadely et al. (2010) 
and Siana et al. (2008) studies, respectively. A2667a and the 
8 O’clock arc have no L(TIR) estimates from far-IR SED fitting 
because the longest band we presently have is 70 /xm. The 70 /xm 
imaging for the latter is complicated by a nearby bright dust- 
obscured galaxy (Section 2.6) that blends with the arc. 

The monochromatic (i.e., single-band) 24 /xm L(TIR) estima- 
tor from Rujopakarn et al. (2012) yields L(TIR) consistent with 
the far-IR SED fitting at 0 < 7 < 2.8 within 0.02 dex, on aver- 
age, with a scatter of 0. 12 dex for individual star-forming galax- 
ies. Rujopakarn et al. (2012) take into account the evolution of 
the bolometric correction from local galaxies out to z = 2.8 (the 
farthest redshift where the 24 /xm band traces predominantly 
dust and PAH emissions) due to the structural differences of 
IR galaxies by using the IR surface brightness, Zzxtir), as an in- 
dicator of IR SED, which has been demonstrated by Rujopakarn 
et al. (201 1). Specifically, high-z IR galaxies are physically ex- 
tended and typically have lOO-lOOOx more surface area than 
local counterparts at similar L(TIR) (Rujopakarn et al. 2011), 
and thus larger PAH-emitting area, which affects the bolometric 
corrections in the PAH wavelength region that is probed by the 
24 /xm band, particularly at z > 1 (Rujopakarn et al. 2012). We 
report the values of L(TIR) from this 24 /xm indicator along 
with those from SED fitting in Table 5. The agreement between 
the monochromatic indicator and the far-IR SED fitting is within 
0.12 dex (average difference of 0.06 dex), similar to the agree- 
ment reported in Rujopakarn et al. (2012). No systematic trends 
in redshift, luminosity, or metallicity are observed. 

The last L(TIR) estimator is derived via the luminosity of 
the 6.2 /xm PAH emission line, following the relationship given 
by Pope et al. (2008), which was measured by combining a 
local sample from Brandi et al. (2006) and the results for their 
submillimeter galaxies (SMGs). Pope et al. (2008) report a slope 
of the relationship between log 6.2 /xm PAH luminosity and log 
L(TIR) to be approximately unity (formal slope = 1.01 ± 0.01) 
over L(TIR) ranging from 10 10 to 10 13 L Q . Our objects are 
shown in Figure 6 on top of data points from Pope et al. 
(2008) and Brandi et al. (2006). The 24 /xm L(TIR) values 
are used here instead of those from the far-IR SED fitting 
because A2667a has no far-IR L(TIR) estimate. Because of 
the methodological differences in the PAH flux measurements 
between Pope et al. (2008) and ours, particularly the assumption 
of the continuum level, for the purpose of the comparison 
in Figure 6, we normalize the PAH luminosity by a factor 
determined by matching our 6.2 /xm PAH luminosity of cB58 
to the measurement of the same object by Siana et al. (2008), 
who adopt the same continuum assumption as Pope et al. (2008). 
The 6.2 /xm PAH feature was chosen to represent the aromatic 
emissions over other stronger PAH emission lines and emission 
complexes because it is narrow and relatively isolated from other 
lines, which mitigates the systematic uncertainties introduced 
by the process of deblending contributions from adjacent lines 
(which is required for, e.g., the 7.7 /xm complex and 8.6 /xm 
line). To estimate the SFR from the 6.2 /xm PAH luminosity, the 
PAH luminosity is first converted to L(TIR) via the Pope et al. 
(2008) conversion and then to SFR via the Rujopakarn et al. 
(2012) relation. 

In addition to estimating SFR via L(TIR), we also show the 
Kennicutt (1998) formula for Ha (no extinction correction) 
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Table 5 

Derived Quantities 


Rujopakarn et al. 


Source 

LIR (log L 0 ) 



SFR (M 0 yr~ 

*) 


a Ha/Spitzer 
/i v 

(8) 

aIR 

A v 

(9) 

A„ 

(10) 

24 /im 
0) 

Far-IR SED 
(2) 

Ha 

(3) 

Ha + 24 /Am 
(4) 

PAH 6 .2 /im 
(5) 

LIR24 /im 
(6) 

LIRpar-IRSED 

(7) 

A2218b 

11.69 

11.68 

10.3 

36.6 

66.7 

48.0 

46.8 

2.8 ± 0.3 

2.1 

2.5 

A2667a a 

11.68 


13.9 

39.9 

19.4 

47.1 


-2.4 a ±3.1 

1.7 

-0.4 

The Clone* 5 

11.20 

11.38 

11.4 

20.0 

36.2 

15.1 

23.2 

<-l.l b 

1.0 

1.0 

A2218a 

11.81 

11.83 

7.5 

42.3 

170.2 

63.8 

66.5 

3.3 ± 0.4 

3.0 

3.2 

A1835a 

12.71 

12.83 

14.5 

298 

558.2 

561 

752 

5.9 ± 0.4 

5.4 

5.7 

cB58 b 

11.04 

11.04 

13.6 

19.8 

26.6 

10.4 

10.4 

<0.7 b 

-0.4 

-0.4 

8 O’clock 

12.17 


71.5 

154 


149 


2.8 ± 0.3 

1.0 

1.9 


Notes. Quantities in this table are corrected for lensing magnification. Column 1: L(TIR) estimated from the monochromatic 24 /Am indicator from 
Rujopakarn et al. (2012). Column 2: L(TIR) measured by fitting SED to far-IR and submillimeter photometry. The values for A2218b and A1835a are 
from Rigby et al. (2008). The values for A2218a, the Clone, and cB58 are from Finkelstein et al. (201 1), Fadely et al. (2010), and Siana et al. (2008), 
respectively. Column 3: Ha SFR estimated by the Kennicutt (1998) formula (no extinction correction). Column 4: extinction-corrected SFR estimated 
using the Ha and rest-frame 24 fim luminosity formula provided by Kennicutt et al. (2009). Column 5: SFR estimated by the PAH 6.2 / im complex 
luminosity. Column 6: SFR estimated by 24 /Am monochromatic indicator provided by Rujopakarn et al. (2012). Column 7: SFR estimated using the 
far-IR L(TIR) in Column 2 and the relationship given by Rujopakarn et al. (2012). Column 8: screen extinction A v in mag estimated with Ha and 
IR recombination lines assuming case-B recombination scenario. The value for cB58 uses Ha and the Paa line, the Bra limit for cB58 indicates an 
extinction limit of A v < 1.8 mag. Column 9: A® in mag. See Section 3.4 for more details. Column 10: final extinction estimate from averaging values 
in Columns 8 and 9. For cB 58 and the Clone, the A^ values are adopted. 
a IR recombination line detected at low significance. 
b IR recombination line not detected. 


10 11 


10 10 


© 


10 8 


10 7 


Figure 6. Relationship of 6.2 /Am PAH emission line luminosity to L(TIR) in 
our sample (large red squares) compared to the local star-forming galaxies 
and LIRGs from the Brandi et al. (2006) study (squares) and the high-redshift 
submillimeter galaxies from Pope et al. (2008, crosses). The dashed line is 
a linear relationship between log luminosities with slope of 1.01 measured 
by Pope et al. (2008). Our galaxies fall within this relationship, confirming the 
approximately unity slope reported by Pope et al. at the intermediate luminosities 
between the local and SMG samples. 

(A color version of this figure is available in the online journal.) 
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Table 6 

Metallicity Estimates from [N n]/Ha 


Source 

Zppo4 

(1) 

Z P P04, 3 rd order 
(2) 

A2218b 

8.55 ± 0.04 

8.53 ± 0.10 

A2667a 

8.67 ± 0.05 

8.73 ± 0.10 

The Clone 

8.49 ± 0.02 

8.45 ± 0.07 

A2218a 

8.55 ± 0.10 

8.52 ± 0.25 

A1835a 

8.75 ± 0.15 

8.91 ± 0.33 

cB58 

8.31 ± 0.10 

8.26 ± 0.35 

8 O’clock 

8.58 ± 0.04 

8.57 ± 0.09 


Notes. Values are in 12+ log(0/H). Column 1: metallicity estimates 
from the N2 linear formula provided by Pettini & Pagel (2004). 
Column 2: metallicity estimates from the N2 third-order estimator 
from Pettini & Pagel (2004). 


can be tested with IR recombination lines that are more robust to 
high extinction (Alonso-Herrero et al. 2006). To obtain a fiducial 
estimate, we assume a foreground screen and measure the 
extinction relative to Ha using the Paa and Bra lines assuming 
case-B recombination. For instance, for Bra, the Ay°^ Spitzei is 
given by 


I Bra ^Bra,0 

= exp 

I Ha ^Ha,0 


. H a/Spitzer 

A v 


1.086 



(1) 


and the Kennicutt et al. (2009) indicator that combines Ha 
and the rest-frame 24 /rm luminosity in order to correct for 
dust extinction in the SFR estimates. The SFR estimates from 
each method are tabulated in Table 5. We discuss the effects of 
extinction on the SFR estimates in Section 4. 

3.4. Extinction 

Extinctions derived from optical lines (e.g., Ha/H/l) could 
suffer systematic underestimation where the overall extinction 
is large and the dust is mixed with the sources. This possibility 


where / Bra ,o/4fo.o = 0.0291 and for Paa, / Paa .,o//Ha,o = 0.1226 
(Osterbrock & Ferland 2006) for T e = 10,000 K at the low- 
density limit. The A Hq: , A Pao , , and A Bra are given by interpolation 
using the Rieke & Lebofsky (1985) extinction law, which for 
A v of 1 mag are 0.8 mag, 0.15 mag, and 0.04 mag, respectively. 

Alternatively, the extinction can be estimated from the optical 
and IR SFR values by taking the latter to be a fiducial SFR. The 
assumption is secure in our sample of IR luminous galaxies, 
where direct UV leakage from the galaxy is expected to be 
small (Rieke et al. 2009, and references therein). This alternative 
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extinction estimate, A ® , following Choi et al. (2006) is 


A 


IR 

Ha 


2.51og 


/ sfr ir \ 

VSFRhJ 


( 2 ) 


The A® can then be converted to the A® via the Rieke & 
Lebofsky (1985) extinction law. When available, we use the 
far-IR SFRir to calculate the A®, otherwise the value estimated 
from the 24 fim L(TIR) is used. 

The J 6^ cl l s P ltze ' method assumes a foreground screen of dust 
and therefore is a lower limit, whereas the A® does not. The 
comparison of measurements from both methods, tabulated 
in Table 5, shows that in three out of four galaxies where 
IR recombination lines are well detected (A2218b, A2218a, 
and A 1835a), both methods agree within the range of uncer- 
tainty, which suggests that the nature of the dust distribution in 
these galaxies is roughly uniform, resembling the foreground 
screen assumption. The other object with a well-detected Paa 
(the 8 O’clock arc), however, shows a 1.8 mag difference be- 
tween a Ho'/Wft''' an£ j ^ir j atter disagreement indicates an 
inhomogeneous mixture of dust in the 8 O’clock arc and high- 
lights the diversity of the dust distribution scenarios at redshift 
1 < z < 3. 

A2218a, A1835a, cB58, the Clone, and the 8 O’clockarc have 
extinction measurements from optical and/or rest-frame optical 
spectroscopy in the literature that can be compared with our 
IR measurements. Richard et al. (2011) found E(B — V) st ar of 
0.18 for A2218a from SED fitting, implying A„ of 0.6 mag, as- 
suming the R = 3.1 law. A1835a has extinction measurements 
by Nesvadba et al. (2007) using Ha/H/J ratio of E(B — V) = 
1.3-1. 6, implying A v of 4. 0-5.0 mag. Teplitz et al. (2000) 
measured E(B — V) of 0.27 for cB58 from the Ha/H/1 
line ratio, implying A v of 0.4 mag. For the Clone, Hainline 
et al. (2009) measured E(B — V ) using the Ha/H y line 
ratio to be 0.28, which implies A v of 0.9 mag. The ex- 
tinction of the 8 O’clock arc was measured by Finkelstein 
et al. (2009) using a weighted mean of Ha, H/3, and H y 
line ratios to be A v = 1.17 ± 0.36 mag (using only 
Ha/H/3 yields E(B — y) gas of 0.97, implying A„ of 1.3 mag); 
and by Dessauges-Zavadsky et al. (201 1) to be E(B — V) gas = 
0.30 ± 0.04 mag, implying A„ of 0.9 mag. From Table 5, the ex- 
tinction measurements from IR recombination lines for A22 1 8a, 
A1835a, and the 8 O’clock arc, where we have secure IRline 
detections, are 3.3 ± 0.4 mag, 5.9 ± 0.4 mag, and 2.8 ± 0.3 mag, 
respectively. While our sample size is too small to draw a general 
conclusion, the comparison suggests that the optical measure- 
ments may sample systematically lower extinction regions, in 
agreement with local LIRGs and ULIRGs (e.g., Alonso-Herrero 
et al. 2006). 

The A ^ a/Sp ' L values for A2667a and the Clone are negative 
because the formal fits to the Paa and Bra line yield values 
(or flux limits) lower than those expected for the extinction of 
Omag (as shown in the simulated A„ in Figures 3 and 4); the 
line is undetected), albeit with a considerable uncertainty. For 
cB58, the line flux for both Paa and Bra are negative and thus 
the A ^°‘/ Sp " ze ' cannot be formally calculated and we report ala 
upper limit as cB58’s j^ cl l s P l,zel _ Given the general agreement 
between j\^ a ' s i n,zel anc | 4 | R j n a majority of objects with secure 
measurements of Paa and Bra, we adopt an average value of 
^Ha/Spii^n as a re p resen t a ti ve value of A„ (except for 

the cB58 and the Clone for which we adopt the A® value). 


CP 

O 

£ 

< 


Figure 7. Optical extinction, A v , from this work compared to the distribution of 
Choi et al. (2006) from the Spitzer First-Look Survey (FLS) with mean redshift 
of z = 0.8. Our points are shown in filled squares, except for the cB58 point 
(empty square), where we adopt the A^ value as an extinction estimate (see 
Section 3.4). The Choi et al. (2006) points are A^ measured by comparing the 
IR SFR, which is taken as a fiducial value, with the optical SFR measured from 
the Ha, H/3, and [On] emission lines. These are shown in red, blue, and gray 
circles, respectively. The intrinsic spread of points toward higher extinction 
of the Choi et al. sample may be larger if the line non-detections (upward 
arrows), which are likely due to large extinction, could have been measured. 
The intrinsic spread of extinction at high-z is far wider than previously known 
from optical-based measurements. 

(A color version of this figure is available in the online journal.) 
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We present our final A„ estimates as a function of L(TIR) 
from the 24 fim indicator in Figure 7. Again, the 24 ji m L(TIR) 
is used here because A2667a and the 8 O’clock arc lack far- 
IR L(TIR) estimates. Our measurements are compared with 
the A® measurements by Choi et al. (2006), who study the 
extinctions of objects selected in the near-IR and mid-IR from 
the Spitzer First Look Survey (FLS), which have a mean redshift 
of z — 0.8 and luminosities ranging from the sub-LIRG to 
ULIRG range. Their extinction measurements are obtained by 
comparing the SFR estimates from the Ha, H p, and [Oil] line 
fluxes with those from IR luminosities. That is, they estimated 
the extinction using the ratio of SFR from the individual optical 
indicators to the IR SFR estimates, e.g., SFRHa/SFRiR, SFRh/i / 
SFRir, etc. The Choi et al. distribution of extinction values 
shows a correlation between A® and L(TIR) with a formal fit 
of A® = 0.75 xlogL(TIR)— 6.35, shown as the orange line in 
Figure 7. The scatter is ~0.8 mag; however, this value does not 
include cases where [O ii] is undetected (shown as lower limits 
in Figure 7). Allowing for these limits, the scatter is >1 mag. 
This fit and the large scatter are similar to the result of Afonso 
et al. (2003), who estimate extinction by comparing Ha and 
[O ii] SFR estimates to those from the 1 .4 GHz radio continuum 
at a similar mean redshift of z = 0.8. A number of authors (e.g., 
Bai et al. 2007; Kocevski et al. 2011) find that the variations 
in extinction are so large that [O ii] can be difficult to detect in 
some galaxies that have high SFRs indicated at 24 /rm, although 
with high-quality spectra there are relatively few cases where 
[O ii] is completely absent (Weiner et al. 2007). 

We found the spread of our extinction values to be at least 
as large. A2218a, A2218b, and the 8 O’clock arc lie within 
the uncertainties of the Choi et al. relationship. A2667a and 
cB58, with nearly zero mag extinction, are less obscured for 


11 


The Astrophysical Journal, 755:168 (13pp), 2012 August 20 


Rujopakarn et al. 




Figure 8. Comparison of star formation rate (SFR) estimated from the uncorrected Ha, Ha + rest-frame L( 24 /im), L(TIR) from 24 /im indicator, and L(TIR) from 
the 6.2 fim PAH emission line, to the SFR estimated from far-IR SED fitted L(TIR), which is taken as a fiducial value, as a function of L(TIR) and A v , shown in the 
left and right panels, respectively. The Ha SFR indicator systematically underestimates L(TIR) with increasing discrepancy from the fiducial at larger L(TIR), while 
the rest are consistent to the fiducial within the range of uncertainties. A trend line for the extinction effect on the Ha SFR indicator in objects with a screen-like 
extinction scenario (Section 4) is shown as the blue dotted line. The points representing the 8 O’clock arc are circled in gray to indicate that L(TIR) of the arc is 
estimated from 24 /im because the SED-fitted L(TIR) is not available. The deviation of Ha SFR of the arc from the trend line in the left panel, however, is not due to 
the L(TIR) measurement, but rather the differences of dust distribution scenarios between the inhomogeneous mixture in the arc and the homogeneous distribution in 
the rest of the sample (see Section 3.4). 

(A color version of this figure is available in the online journal.) 


their L(TIR), whereas the extinction of A1835a is far above the 
relationship. Selection effects may have reduced the scatter in 
many previous studies (e.g., the need to have an [Oil] line of 
measurable strength). Our results indicate that galaxies up to 
z ~ 3 show a very large range of A v , as is the case in the local 
universe. 

4. EFFECTS OF EXTINCTION ON STAR FORMATION 
RATE INDICATORS 

In Figure 8, we compare the SFR estimates from Ha (un- 
corrected for dust extinction), Ha + L( 24 /rm), L(24 //in), and 
the 6.2 /rm PAH emission with those of L(TIR) measured by 
far-IR SED fitting as a function of L(TIR) (Figure 8, left) and 
A v (Figure 8, right). A2667a is excluded from this comparison 
due to its lack of the far-IR SED fitted L(TIR) and the fact that 
the curvature of the lensed arc caused a fraction of Ha flux to 
fall off the slit (the position of the slit is shown in the bottom 
left panel of Figure 1; a part of the “A2” arc is outside the slit). 
Without the correction for dust extinction, the Ha SFR indicator 
is clearly affected at high L(TIR) and extinction. For the cases 
that extinction behavior resembles a uniform dust screen (every 
galaxy except the 8 O’clock arc), there is a tight trend of the Ha 
SFR deviation from the fiducial SED-fitted SFR as a function 
of the log[L(TIR)] that is given by 

log[SFR(Ha)/SFR(LIR)] = 11.21 - 1.011og[L(TIR)]. (3) 

Applying a single overall extinction to Ha will increase the 
estimates of SFRs and improve the agreement to the L(TIR) 
SFR, but will not correct the trend with L(TIR). Further 
improvement in Ha-based SFRs require introducing corrections 
for extinction as a function of stellar mass (e.g., Moustakas et al. 
2006; Weiner et al. 2007), or of the SFR itself (e.g., Buat et al. 
2005); see Section 4. 

The deviation of the Ha point of 8 O’clock arc in 
Figure 8 (left) indicates that a larger amount of Ha flux is 


escaping from the arc given its L(TIR), which could be due 
to a different dust distribution relative to the star-forming re- 
gions. The indicator is usually applied with a single nom- 
inal level of extinction; in this case, the line in Figure 8 
(left) would be shifted upward but the slope would not be 
corrected. 

For A2218a, A 1835a, and the 8 O’clock arc where Paa line is 
well detected, we apply the extinction correction to the Paa line 
luminosity and compare the ratio of L(24/xm)/L(Paa corr ) with 
the relationship that Alonso-Herrero et al. (2006) found for local 
LIRGs. Papovich et al. (2009) made this comparison for A2218a 
and found the L(24/xm)/L(Paa corr ) ratio to be ~ 0.5 dex 
lower than the local relation, indicating that A2218a has lower 
L(24 |Um) than local galaxies of comparable L(Paa corr ), while 
agreeing with those of local individual H ii regions. Papovich 
et al. interpreted this difference as an indication that A2218a 
harbors extended star-forming regions similar to a scaled- 
up local Hii regions, rather than nuclear starburst like local 
U/LIRGs. In our analysis, we have reproduced the IA24 //in)/ 
L(Paa corr ) ratio measured by Papovich et al. using our reduction 
technique for A2218a, log[L(24 /im)/L(Paa corr )] = 2.2 ± 0.2, 
and also found that the ratio for A1835a and the 8 O’clock arc 
are 0.9-1. 1 dex lower than the Alonso-Herrero et al. (2006) 
relationship. This is consistent with the result in Table 5 that the 
Rujopakarn et al. (2012) single-band 24 fim L(TIR) indicator 
yields L(TIR) values that agree with the far-IR SED-fitted 
L(TIR), which implies that these galaxies have extended star 
formation (see Section 3.3 for details of the indicator). That is, 
the local relationship of L(Paa corr ) and L( 24 /rm) reported by 
Alonso-Herrero et al. (2006) could have limited applicability at 
high redshift because the Paa line will systematically have lower 
extinction for a given L(24 // m ) due to the extended structure 
of star formation. 

All the other indicators give consistent estimates, within the 
expected errors of ~0.2 dex. We found that the Ha + L(24 /imj 
indicator by Kennicutt et al. (2009) tends toward a smaller value 
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of SFR at large extinctions or large L(TIR), particularly above 
10 12 L q , as observed in A1835a and the 8 O’clock arc. The 
SFR estimates from the 6.2 /rm PAH emission do not show a 
systematic trend with L(TIR), although the scatter is larger 
than with other indicators. The scatter for individual galaxies 
is ~ 0.2 dex, similar to the scatter of L(TIR) values we have 
found in Rujopakarn et al. (2012) at redshifts where the 24 /im 
band probes the PAH emissions. Since these are high S/N 
measurements of an individual (i.e., relatively non-blended) 
PAH emission line, this result suggests that the Rujopakarn et al. 
(2012) SFR indicator has succeeded in correcting for the SED 
evolution. 

5. CONCLUSIONS 

We observed four strongly gravitationally lensed star-forming 
galaxies at 1 < z < 3 with Spitzer/IRS and the LBT / 
LUCIFER to obtain the mid-IR and near-IR spectroscopy. 
These observations are targeted to cover IR recombination lines, 
including Ha in the near-IR and Paa or Bra in the mid-IR. We 
include another three galaxies from the literature with similar 
suites of observations, yielding a total sample of seven galaxies. 
Our sample spans the redshift range of 1 . 03-2 . 7 3 and the L (TIR) 
range of 1.3 x 10 n L Q to 7.0 x 10 12 L Q . 

The IR recombination line ratios are used to measure extinc- 
tion that can probe deep into the highly obscured star-forming 
regions and thus provide an unbiased measure of extinction 
under the foreground screen assumption. Independently, we es- 
timate the extinction by comparing the optical and IR SFRs, 
a method that does not make the foreground screen assump- 
tion. The results from the two methods are consistent in three 
out of four galaxies with good IR recombination line flux mea- 
surements, suggesting that the dust extinction in these galax- 
ies is consistent with a foreground screen (i.e., uniform dust 
distribution). However, in the fourth case, the extinction es- 
timates from two methods disagree by 1.8 mag, indicating a 
deviation from the uniform dust screen assumption, which sug- 
gests an inhomogeneous dust mixture. The extinction range of 
our sample (assuming a foreground screen) is ~0. 0-5. 9 mag, 
which is a larger spread than previously known for intermediate- 
and high-redshift galaxies based on measurements with opti- 
cal emission lines. These results suggest a large diversity in 
both the extinction levels and dust distribution scenarios at high 
redshift. 

We compare the performance of various SFR indicators 
over the extinction range and find that substantial extinc- 
tion corrections are required for the Ha-based SFR indica- 
tor. The remaining indicators (1) combined Ha and L(24/rm), 
(2) L( 24 /rm), and (3) PAH (6.2 /im) all give estimates consis- 
tent to within the expected uncertainties of ~0.2 dex. 

We thank Brian Siana and Eiichi Egami for helpful discus- 
sions. W.R. thanks Alexandra Pope and Philip Choi for data 
points in Figures 6 and 7, respectively, and acknowledges the 
support from the Thai Government Scholarship. This work is 
supported by contract 1255094 from Caltech/JPL to the Uni- 
versity of Arizona. 


REFERENCES 

Afonso, J., Hopkins, A., Mobasher, B., & Almeida, C. 2003, ApJ, 597, 269 
Allam, S. S., Tucker, D. L., Lin, H., et al. 2007, ApJ, 662, L51 
Alonso-Herrero, A., Rieke, G. H., Rieke, M. J., et al. 2006, ApJ, 650, 835 
Bai, L., Marcillac, D., Rieke, G. H., et al. 2007, ApJ, 664, 181 
Becker, G. D., Rauch, M., & Sargent, W. L. W. 2009, ApJ, 698, 1010 
Berta, S., Magnelli, B., Nordon, R., et al. 2011, A&A, 532, A49 
Bian, F., Fan, X., Bechtold, J., et al. 2010, ApJ, 725, 1877 
Bouche, N., Lehnert, M. D., Aguirre, A., Peroux, C., & Bergeron, J. 
2007, MNRAS, 378, 525 

Brandi, B. R., Bernard-Salas, J., Spoon, H. W. W., et al. 2006, ApJ, 653, 1 129 
Buat, V., Iglesias-Paramo, J., Seibert, M., et al. 2005, ApJ, 619, L51 
Choi, P I., Yan, L., Im, M., et al. 2006, ApJ, 637, 227 

Dessauges-Zavadsky, M., Christensen, L., D’Odorico, S., Schaerer, D., & 
Richard, J. 2011, A&A, 533, A15 
Diamond-Stanic, A. M., & Rieke, G. H. 2010, ApJ, 724, 140 
Diolaiti, E., Bendinelli, O., Bonaccini, D., et al. 2000, A&AS, 147, 335 
Dole, H., Lagache, G., Puget, J.-L., et al. 2006, A&A, 451, 417 
Engelbracht, C. W., Blaylock, M., Su, K. Y. L., et al. 2007, PASP, 1 19, 994 
Erb, D. K., Shapley, A. E., Pettini, M., et al. 2006, ApJ, 644, 813 
Fadely, R., Allam, S. S., Baker, A. J., et al. 2010, ApJ, 723, 729 
Finkelstein, K. D., Papovich, C., Finkelstein, S. L., et al. 201 1, ApJ, 742, 108 
Finkelstein, S. L., Papovich, C., Rudnick, G., et al. 2009, ApJ, 700, 376 
Hainline, K. N., Shapley, A. E., Kornei, K. A., et al. 2009, ApJ, 701, 52 
Kennicutt, R. C. 1998, ARA&A, 36, 189 

Kennicutt, R. C., Jr., Hao, C.-N., Calzetti, D., et al. 2009, ApJ, 703, 1672 
Kneib, J.-R, Ellis, R. S., Smail, I., Couch, W. J., & Sharpies, R. M. 1996, ApJ, 
471,643 

Kneib, J.-R, van der Werf, P. R, Kraiberg Knudsen, K., et al. 2004, MNRAS, 
349,1211 

Kobayashi, C., Springel, V., & White, S. D. M. 2007, MNRAS, 376, 1465 
Kocevski, D. D., Lemaux, B. C., Lubin, L. M., et al. 2011, ApJ, 736, 38 
Kroupa, P. 2002, Science, 295, 82 

Le Floc’h, E., Papovich, C., Dole, H., et al. 2005, ApJ, 632, 169 
Lin, H., Buckley-Geer, E., Allam, S. S., et al. 2009, ApJ, 699, 1242 
Magnelli, B., Elbaz, D., Chary, R. R., et al. 2011, A&A, 528, A35 
Mandel, H., Seifert, W., Lenzen, R., et al. 2007, Astron. Nachr., 328, 626 
Moustakas, J., Kennicutt, R. C., Jr., & Tremonti, C. A. 2006, ApJ, 642, 775 
Nesvadba, N. P. H., Lehnert, M. D., Genzel, R., et al. 2007, ApJ, 657, 725 
Osterbrock, D. E., & Ferland, G. J. 2006, in Astrophysics of Gaseous Nebulae 
and Active Galactic Nuclei, ed. D. E. Osterbrock & G. J. Ferland (2nd ed.; 
Sausalito, CA: Univ. Science Books) 

Papovich, C., Rudnick, G., Rigby, J. R., et al. 2009, ApJ, 704, 1506 
Perez-Gonzalez, P. G., Rieke, G. H., Egami, E., et al. 2005, ApJ, 630, 82 
Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 
Pope, A., Chary, R.-R., Alexander, D. M., et al. 2008, ApJ, 675, 1171 
Reddy, N. A., Steidel, C. C., Pettini, M., et al. 2008, ApJS, 175, 48 
Richard, J., Jones, T., Ellis, R., et al. 2011, MNRAS, 413, 643 
Rieke, G. H., Alonso-Herrero, A., Weiner, B. J., et al. 2009, ApJ, 692, 556 
Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618 
Rieke, G. H., Young, E. T., Engelbracht, C. W., et al. 2004, ApJS, 154, 25 
Rigby, J. R., Marcillac, D., Egami, E., et al. 2008, ApJ, 675, 262 
Rodighiero, G., Vaccari, M., Franceschini, A., et al. 2010, A&A, 515, A8 
Rujopakarn, W., Eisenstein, D. J., Rieke, G. H., et al. 2010, ApJ, 718, 1 171 
Rujopakarn, W., Rieke, G. H., Eisenstein, D. J., & Juneau, S. 2011, ApJ, 726, 
93 

Rujopakarn, W., Rieke, G. H., Weiner, B. J., et al. 2012, arXiv:l 107.2921 
Sand, D. J., Treu, T., Ellis, R. S., & Smith, G. P. 2005, ApJ, 627, 32 
Seitz, S., Saglia, R. P, Bender, R., et al. 1998, MNRAS, 298, 945 
Siana, B., Teplitz, H. I., Chary, R.-R., Colbert, J., & Frayer, D. T. 2008, ApJ, 
689, 59 

Smail, I., Smith, G. P, & Ivison, R. J. 2005, ApJ, 631, 121 
Smith, J. D. T., Draine, B. T., Dale, D. A., et al. 2007, ApJ, 656, 770 
Teplitz, H. I., Desai, V., Armus, L., et al. 2007, ApJ, 659, 941 
Teplitz, H. I., McLean, I. S., Becklin, E. E., et al. 2000, ApJ, 533, L65 
Vacca, W. D., Cushing, M. C., & Rayner, J. T. 2003, PASP, 1 15, 389 
Weiner, B. J., Papovich, C., Bundy, K., et al. 2007, ApJ, 660, L39 
Yee, H. K. C., Ellingson, E., Bechtold, J., Carlberg, R. G., & Cuillandre, J.-C. 
1996, AJ, 111, 1783 


13 


